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Abstract. We present a renormalization approach to solve the Sznajd opinion formation model on complex networks. 
For the case of two opinions, we present an expression of the probability of reaching consensus for a given opinion 
as a function of the initial fraction of agents with that opinion. The calculations reproduce the sharp transition of the 
model on a fixed network, as well as the recently observed smooth function for the model when simulated on a growing 
complex networks. 
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1 Introduction 

Opinions can either be made up by a person or taken over from 
another person. Sometimes some people try to force their opi- 
nions on others. In general, all people are free to form opinions 
as they see fit. The mechanism of opinion formation is "norma- 
tive", i.e., normative in the sense of what ought to be, opposed 
to a "positive" mechanism, which is based on observation what 
is 1 1 1. Based on this facts, and with the necessary simplifying 
assumptions, socio-physics gave the opportunity to apply tech- 
niques of statistical physics to model opinion formation among 
people 1 2I3EI . 

One of the opinion formation models that has generated 
immediate interest in many authors on the field is the Sznajd 
model 121, which is based on the slogan "together we stand": 
Individuals are represented by the lattice nodes (one-dimensional 
in its first version), and each randomly selected pair of neigh- 
bors convinces all their neighbors of their opinions, if and only 
if the pair shares the same opinion; otherwise, the neighbors' 
opinion are not affected. It differs from other consensus mod- 
els by dealing only with communication between neighbors, 
and the information flows outward as in rumor spreading: a 
site does not follow what the neighbours tell the site. 

On networks with fixed size, the results of the model do not 
depend much on the spatial dimensionality and type of neigh- 
borhood selected (i.e., two nodes convince the others, three 
convince the others, etc.) |6 7 9|. In the case of q choices of 
opinion, the model has q homogeneous absorbing states, where 
all individuals choose the same opinion; in the context of opin- 
ion, one says the system reaches consensus. The case of two 
opinions {q ~ 2) has been the most studied, denoting opinions 
as Ising variables "up" or +1, and "down" or —1. In more than 



one dimension, the probability {Pup) of reaching consensus "all 
up" depends on the initial fraction p of individuals with opinion 
"up"; for p > 0.5, the probability of reaching "all up" as sta- 
tionary state is close to one, while forp < 0.5 it is negligible, 
having a sharp transition inp = 0.5, which can be interpreted 
as a dynamical phase transition. Computer simulations in |8| 
indicate that the universality class associated with this dynam- 
ical phase transition is different from the universality class of 
the Ising model. The distribution of time needed to reach the 
stationary state is a peak followed by a fast decay flO]. 

Much less is known about the Sznajd model on growing 
networks. Interactions of groups of people in some circum- 
stances can be thought as a growing system, i.e., in a city with 
positive rate of immigration. In a first and simple approxima- 
tion, it can be modeled by a growing scale-free network II II . 
Recently, applying a Sznajd model recipe not after the com- 
plete network has been constructed, but while the network grows, 
i.e, while each new node is added to the network, one could ob- 
serve that the Sznajd model simulated on scale-free networks, 
Barabasi- Albert network and a pseudo-fractal network 1121 . the 
system reaches consensus fTTT4l. But in contrast to the sharp 
transition observed for the networks of fixed size, in which the 
Sznajd recipe is performed only after building up completely 
the network, the probability that the system reaches "all up" 
for a growing complex networks is a smooth function of p. In 
addition, this function depends on the type of neighborhood 
selected. 

In this work, we propose a real space renormalization ap- 
proach 1 15 1 to calculate the probability P^p (p) of reaching con- 
sensus on opinion "up" as a function of the initial fraction 
p of opinion "up". Our results are for two common rules of 
neighborhood, namely "r-convince all their neighbors", with 
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r ~2 and r = 3. We have obtained the two well-known results 
known for the model: a smooth function of p for the growing 
case and an expression which approximates the step function 
for fixed networks. 

In the next section, we present the hierarchical network 
used in our calculations. Then, we present the renormalization 
approach and the analytical expressions obtained, each case is 
compared with the results from the numerical simulations, pre- 
viously reported in 1141 . as well as for the BA scale-free net- 
work. 




t=0 t=l t=2 

Fig. 1. The first three generations of the scale-free pseudo-fractal 
graph. At each iteration step t, every edge generates an additional ver- 
tex, which is attached to the two vertices of this edge. 



2 Hierarchical Networit 



to a chosen Sznajd rule (r) in a growing network {g). Subse- 
quent self-iterations of Pup{p)\r,g result in a step function, i.e., 
Pup{p)\r.g — Pup{p)r.f, where the subscript index / corre- 
ponds to the result obtained for a network of fixed size. 

For a population fraction p with opinion "up", the general 
method is as follows: 

- Given a neighborhood rule r, the chosen basic cell corre- 
sponds to the minimum generation t of the hierarchical net- 
work, such that r > Nt (the r agents must have at least one 
agent to convince). We call this resulting number of nodes 
in the cell rir- 

- The probability of each possible configuration in a elemen- 
tary cell is easily calculated, such that 

1 = Pall{p)\r = ^S„^fe/(1 -p)"'-^ (2) 
k=2 

with the binomial coefficient Bn^k over the appropriate cell 
for the chosen rule: 

S„,fc = n,!/[fc!K-fc)!]. (3) 

- From all the configurations calculated above, we select the 
subset that gives "all up" when applying the selected Sznajd 
rule on the cell, the sum of all of them is Pup{p)\r.g- 

Pup{p)\r.g ^ Pau{p)\r,up (4) 

Next, we illustrate the result of the method with r = 2 and 
r = 3. 



The deterministic scale-free graph used in this work grows as 
follows: At each time step, every edge generates an additional 
vertex, which is attached to both end vertices of this edge. Ini- 
tially, at t = 0, we have a triangle of edges connecting three 
vertices, at t = 1, the graph consists of 6 vertices connected by 
9 edges, and so on (see FigQ. The total number of vertices at 
iteration t is 

3(3* + 1) 



(1) 



This simple rule produces a complex growing network. Such a 
graph is called a pseudo-fractal. In the next section, we present 
the use of this hierarchical network to find expressions that 
agree with the simulated results of the Sznajd model on com- 
plex networks. 



3 Renormalization Approach 

Our method can be very intuitive and is based on the method 
proposed by Galam to study bottom-up democratic voting by 
majority rule in a square lattice |3|, where the predictions of 
the results in all the lattice are based on the applications of the 
majority rule over a basic cell of neighbors, called renormal- 
ization cell. 

We find that given a neighborhood rule, it is enough to 
choose an appropriate generation of a hierarchical network for 
calculating P^p {p)\r.g, which agrees with the the numerical re- 
sults of the model on growing networks The subscript index 
r, g in Pup {p)\r.g is to stress that the resulting function belongs 



*BA nstwork (growing to N|=29526) 
^ - ■♦Pseudof racial (growing to Nt=29576) 
RSRG calculation (sq. (6)) 




0.4 0.6 
initial density (p) 

Fig. 2. Comparison between the function presented in Eq. |6| (solid 
line) with Monte Carlo simulations on a growing pseudo-fractal (trian- 
gles with error-bars) and on a growing BA scale-free network (stars). 
In both networks, 29576 nodes are considered. We count the number 
of samples, out of 1000, for which the fixed point all "up" is obtained 
when different values for the initial concentration p of nodes "up" are 
simulated for rule r — 2. 
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3.1 Case r = 2 

3.1 .1 Growing 

For r = 2, the triangle of the generation t = is the basic 
cell. Thus n,. = 3 and, for a given fraction p, all the possible 
configurations are: 

1 = Pauip)\2 =p' + 2p\l -p) + 2p{l -pf + [l-pf (5) 

If we apply the selected Sznajd rule r = 2 over the triangle, 
only the configurations expressed in the first two terms of the 
sum give "all up". Therefore: 



P^.p{p)\2,g - + 2p2(l - p) = ip^ _ 2p^ 



(6) 



In Fig.|2j we can see the good agreement of Eq.|6lwith the nu- 
merical results 1 14 1 for the Sznajd model on a growing pseudo- 
fractal, as well as for the Barabasi-Albert scale-free network 

CD. 



3.1.2 Fixed 



^ - .<> Sznajd model (fixed Nt=29576) 
B - d Sznajd model (fixed Nt=7971 63) 
O - €>Sznajd model (fixed Nt=2391486) 
RSRG calculation (n.=100, for eq. (7)) 
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Fig. 3. Eq. 0with rii = 100 (solid line) compared with simulations 
on a fixed pseudo-fractal networks with Nt — 29576, 797163 and 
2391486 nodes (dashed line with symbols). Other simulation condi- 
tions as presented in the caption of Fig.|2| 



In order to recover the reported result on a fixed network, 
one makes renormalization iterations, which means simply self- 
composing the Eq.|6j 



Pup{p)\2J ^ P:;{p)\2,g, 



(7) 



and in the limit of large number of iterations (rii — 1), one re- 
covers the step function observed numerically for the model 



on fixed networks. Note that the number of terms and the co- 
efficients sizes increase very fast, as one can observe in the 
expression of only one composition: 

p2^(p)|2 g = 27 p^ - 36p5 - 42/ + 108/ - 72/ + 16/, 

(8) 

therefore, the multiple compositions presented in Fig.|2]are it- 
erated with a computer. Figure|3]shows that the numerical sim- 
ulations on large networks tend to the step function calculated 
fromEq.Hwithrij = 100. 



3.2 Case r 3 

3.2.1 Growing 



^BA network (growing to N,=29526) 
A - -A Pseudofractal (growing to Nt=29526) 
RSRG calculation (eq. (10)) 




0.4 0.6 
initial density (p) 

Fig. 4. Eq. 1101 (solid line) compared with the results from the sim- 
ulations on a growing pseudo-fractal (triangles with error-bars) and 
on a growing BA scale-free network (stars) when r = 3. The other 
simulation conditions are the same of Fig.|2| 



The core of the method is the selection of the correct con- 
figurations after applying the Sznajd rule on it. As we will see 
for this rule, when the number of nodes in the renormalization 
cell is even, there are some symmetrical configurations which 
can have either "all up" or "all down" with the same proba- 
bility. In this case only half of them are summed to Pup- For 
r = 3, the generation t = 1 is the basic cell. Thus Ur = 6 and, 
for a given fraction p, all the possible configurations are: 



I ^ Pau{p)\3 = (I + {I - P) f 



(9) 



Note that the values of the binomial coefficient in the consecu- 
tive terms are: 1, 6, 15, 20, 15, 6, 1. From the 20 configurations 
of the ith term, there are 7 that give "all up"(shown in Fig.|6l 
at Appendix|3, the coiTesponding 7 opposed cases which give 
"all down", and 6 symmetrical configurations shown in Fig.0 
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(Appendix]^ that can give either "all up" or "all down". There- 
fore, these group of configurations contribute with 7 + 0.5 x 6, 
and we have: 

Pup{p)kg = P^'+Gp' (1 - p) + 15/ (1 - p)' + 10/ (1 - pf 

(10) 

In Fig.|4] we see that Ea.llOlagrees very well with the nu- 
merical results 1 14 1 for the Sznajd model on a growing network 
when the rule r = 3 is considered. 



3.2.2 Fixed 



O - -O Sznajd model (fixed, Nt=29526) 
B - €] Sznajd model (fixed, Nt = 88575) 
^ - -OSznajd model (fixed, Nt.797163) 
RSRG (eq. (10) self-composed witfi n,:=10) 
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Fig. 5. Eg . 1 1 Ol self-composed in 9 steps (solid line) compared with the 
result from the simulations on a fixed pseudo-fractal networks with 
r — 3 (dashed lines with symbols). The other simulation conditions 
are the same as presented in the caption of Fig.|3| 



The result of the composition for this case is far more com- 
plicated and only 1 self-composition of eq.^](ni = 2) akeady 
needs a computer, as shows the following expression: 



Pupipy^ 



3,3 



-1249989 /2_,_390897pii. 



158184pi°+28561p' 



270741222 pi^-lGG735317p^^+41109081 



643783179 p 

- 17504838 + 5585931 - 15244686567^^4+ 

11863411551 p^^ - 7642674243 p^^ + 4315583718 p^^- 

2347570026 + 1281132990 - 816731505^^"+ 

2281401855p29 - 5100164190^^^ + 9199907505 p^^- 

13440029166 p2^ + 15908268375 p^^ - 2187p3^ + 65610 p^^ 

-925101 p^'^ + 8148762 p^^ - 50268195 p^^ + 230706630 p^^ 

In Fig.|5]we see the step function obtained with only 9 steps 
of composition compared with the numerical results on a fixed 
network of different sizes; as we see the results agree very well 
with the simulations of the model on large networks. 



4 Conclusions 

Based on opinion formation rules of the usual Sznajd model, 
we use a renormalization approach to give an expression for 
the probability of consensus into one opinion as a function of 
the initial fraction of this opinion. 

We show that for a given Sznajd rule it is enough to solve 
exactly the model on an appropriate basic cell in order to find 
an expression for the smooth function, found numerically for 
the model on a growing network. Several self-compositions of 
the obtained function give the step function observed for the 
model on a network of fixed size. Further renormalization pat- 
terns has to be tested, but in order to reproduce the results of 
the Sznajd model on growing SF networks, a SF hierarchical 
network must be chosen. 

The proposed method could be, in principle, extended to 
other types of neighborhood and more interestingly to many 
choices of opinion (q > 2), which is an feature of the model 
used to simulate elections processes L16,..14 10 1. obtaining re- 
sults consistent with some empirical observations lfTTl . 
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Appendix A: Configurations witli tlie same fraction of 
nodes "up" and "down" 

Here we present some of the possible configurations apply- 
ing the Sznajd rule corresponding to the r = 3 on its appro- 
priate renormalization pattern (n^ = 6). In particular, we show 
the case of half of the nodes having opinion up, mentioned in 
Section l3.2.1l and represented by the fourth term in Eg. 1101 

Figure |6l shows the 7 configurations that give as a result 
"all up", when applying the Sznajd rule, i.e., three consecutive 
nodes with opinion +1 convince all their neighbors. Note that 
interchanging + and — , we have the 7 configurations for the 
opposed case of consensus "all down". 

Figure0presents the 6 symmetrical configurations that have 
cSconsecutive nodes with +1, as well as 3 nodes with —1 giv- 
ing consensus "all up" or "all down", respectively. Thus, these 
configurations contribute wit 0.5 x 6 to the probability of con- 
sensus "all up"(Pp), as showed in Section l3.2.1l 
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Fig. 7. Configurations that generate either consensus "all up" or "all 
down" with r = 3. 
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